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Abstract 

We formulate higher order variations of a Lagrangian in the geomet- 
ric framework of jet prolongations of fibered manifolds. Our formalism 
applies to Lagrangians which depend on an arbitrary number of indepen- 
dent and dependent variables, together with higher order derivatives. In 
particular, we show that the second variation is equal (up to horizontal 
differentials) to the vertical differential of the Euler-Lagrange morphism 
which turns out to be self-adjoint along solutions of the Euler-Lagrange 
equations. These two objects, respectively, generalize in an invariant way 
the Hessian morphism and the Jacobi morphism (which is then self-adjoint 
along critical sections) of a given Lagrangian to the case of higher order 
Lagrangians. Some examples of classical Lagrangians are provided to il- 
lustrate our method. 
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1 Introduction 

An important aspect of mathematics that can be fit into differential geometry 
is the calculus of variations. This research started with several formulations 
of calculus of variations on jet spaces (see, e.g. , [151 l2~2"] and Appendix 2). 
Jet spaces are the natural framework for differential equations in differential 
geometry. Good sources are g| EI3 EH EID E<3 E3 E3| It become later evident 
that the passage from Lagrangians to Euler-Lagrange equations was nothing but 
a differential of a certain complex [341 1351 1361 1371 EH| : this lead to variational 
sequences and much more. 

So far we found that an interesting aspect of calculus of variations was not 
developed in all details from the point of view of geometric formulations on 
jet spaces: namely, the second and higher variations of Lagrangians. The La- 
grangian characterization of the second variation of a Lagrangian in the frame- 
work of jet bundles has been considered in ^H] and, more recently, in 0E1EIE3- 
We stress that in [Tll%l ll5| only first-order Lagrangians were considered, while in 
[H] also a distinguished class of second-order Lagrangians has been studied. In 
particular, in [3E1EI it was shown how to recast (up to divergencies) the system 
formed by the Euler-Lagrange equations together with the Jacobi equations for 
a given Lagrangian as the Euler Lagrange equations for a deformed Lagrangian. 

In this paper we provide a geometrical characterization of the second vari- 
ation of a Lagrangian of arbitrary order in the general case of n independent 
variables and m unknown functions. The second variation can be written in an 
infinite number of ways, by adding arbitrary total divergencies. Some prelimi- 
nary results concerning the Jacobi morphism were obtained in [111 1121 1^5] , In 
the present paper, we exhibit a distinguished representative in the class of all 
such forms. Our representative has remarkable intrinsic and coordinate inter- 
pretations. In particular, we prove that, within the framework of finite order 
variational sequences, the Jacobi morphism turns out to be self-adjoint along 
critical sections, i.e. , along solutions of the Euler-Lagrange equations (see also 
[Sn] for some important consequences of this fact). We show clearly the connec- 
tion between Hessian and Jacobi morphisms; furthermore, their representatives 
are ready for applications to arbitrary order Lagrangians, as we show in the 
examples. 

Notice that in literature only the first order case is usually treated: see 
e.g. |1)JEI1E1E3- Moreover, our formulation has the advantage to be easily 
generalizable to higher variations. In ^Oj a Poincare-Cartan form is used in the 
first-order case and for one independent variable, but in higher-order calculus of 
variations the Poincare-Cartan form is no longer unique, so that the approach 
would not lead to a unique formulation of variations of higher order Lagrangians. 
Last but not least, our approach can be further generalized to higher variations 
of Euler-Lagrange type morphisms, Hclmholtz morphisms and to all forms of 
the variational sequence. This is relevant also in view of the role played by 
the second variation in many geometric contexts. For example, the second 
variation of the Yang-Mills functional has an algebraic structure which leads 
in dimension 4 to important geometric consequences about stable Yang-Mills 
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connections, such as local minima of the functional H |SJ. As well as, it is 
relevant for the study of the theory of stable and unstable minimal submanifolds 
of a Riemannian manifold [32] . 

The paper is organized as follows. In section 2 we recall elements of calculus 
of variations on jet spaces. We use the language of finite order variational 
sequences, as developed by Krupka |23| . 

In section 3 we introduce the notion of i-variation of a section as an i- 
parameter 'deformation'. The variations that we consider are of general (non- 
linear) form 23]. This allows us to pass to infinitesimal variations, i.e. Lie 
derivatives, in a natural and straightforward way. Then, we introduce variation 
of forms as derivatives of their pull-back through the variation of a section with 
respect to the parameters. 

Such a notion of variation is applied in section 4. We concentrate ourselves 
on the computation of the second variation of a Lagrangian. It turns out that, 
on critical sections (i.e. on solutions of the Euler-Lagrange equations), the sec- 
ond variation equals (up to total divergencies) the vertical differential of the 
Euler-Lagrange morphism as well as its adjoint morphism. Thus we first of all 
generalize to higher order Lagrangians the well known fact (in the first order 
case) that the second variation coincides with the Hessian morphism up to total 
divergences (see, e.g. ^] El); furthermore we characterize the Jacobi mor- 
phism of a given Lagrangian as the vertical differential of the Euler-Lagrange 
morphism, which we show here to be self-adjoint along critical sections. We 
stress that, for a first order Lagrangian, these two morphisms coincide with the 
standard Hessian and Jacobi maps. Hence, our method generalizes maps related 
to the second variation to the case of arbitrary order Lagrangians. 

In the last section we consider a few well known and significant examples, 
showing the role played in these cases by the geometric objects above. 

Two short Appendices are included, in which we provide a synthetic version 
of the well-known jet space formulation of variational problems, both for con- 
venience of the less experienced reader and for a better understanding of our 
motivation and formalism. 

2 Variational sequences on jets of fibered man- 
ifolds 

We recall in this section some basic definitions and results from the theory of jet 
spaces. Complete treatments of this subject, with different characters, can be 
found in d |13 El H2 ESI Our exposition follows more closely HO Eg]. 

Our framework is a fibered manifold ir: Y — ► X, with dimX = n and 
dim Y = n + m. We recall that a fibered manifold is just a surjective submersion 
7r; in other words, fibers of 7r need not to be mutually diffeomorphic. A section 
of 7r is defined to be a map s : X — > Y such that it o s — idx ■ We denote by 
VY :=kerT7r C TY the vertical subbundle of the tangent bundle TY. 
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2.1 Jet spaces 

For r > we are concerned with the r-jet space J r Y . This space is defined as 
the set of equivalence classes of sections of tt having a contact of order at least 
r at a given point. Equivalent sections have the same p-th order differential at 
the given point, p < r. We set J Y = Y. 

There are natural projections ir r s : J r Y — ► J S Y, r > s, sending r-th equiva- 
lent sections into s-th equivalent sections. Moreover, there are obvious natural 
projections tt t : J r Y — > X. The spaces J r Y are endowed with a differentiable 
structure making 7r£ fiber bundles. Among these, it can be proved that tt^._ 1 
are affine fiberings (r > 1). 

Charts on Y adapted to tt are denoted by (x x ,y i ). Greek indices A, (i, . . . 
run from 1 to n and they label base coordinates, while Latin indices . . . run 
from 1 to m and label fiber coordinates, unless otherwise specified. We denote 
by (d\, di) and (d x , d z ) the local bases of vector fields and 1 -forms on Y induced 
by an adapted chart, respectively. 

We denote multi-indices of dimension n by boldface Greek letters such as 
er = (di, . . . ,cr n ), with < cr p , /i = l,...,ri; we set \<r\ :=<7i + ■ • • + cr„ and 
ct\:=o\ \ . . .o n \. The charts induced on J r Y are denoted by {x x ,y l rT ), with < 
\<r\ < r; in particular, we set y % Q = y % . The local vector fields and forms of J r Y 
induced by the above coordinates are denoted by (df ) and (d^), respectively. 

We recall that a section s : X — > Y can be prolonged to a section j r s : X — > 
J r Y. If we set y % o s = s 4 , then we have the coordinate expression 

If the order of prolongation needed for s in formulae is clear from context, then 
we simply denote the prolongation of s by js. 

The jet spaces carry a natural structure, the Cartan ( or contact) distribution 
4 . It is the vector subbundle of TJ r Y generated by vectors which are tangent 
to submanifolds of the form j r s(X) C J r Y. We present here a variant of this 
structure (see e.g. |25|). 

We consider the natural complementary fibered morphisms over the affine 
fibering J r+ \Y — > J r Y induced by contact maps on jet spaces 

D : J r+1 Y x TX -> TJ r Y, u : J r+1 Y x TJ r Y VJ r Y, 

X J r Y 

with coordinate expressions, for < \<r\ < r, given by 

D =d x ®D x = d x ®(d x + y j <T+x dj), lo =wi®dj = (dl-yi +x d x )®dj. 

Here, the map D is the inclusion of TX into TJ r Y through the differential Tj r s 
of any prolonged section j r s, while lo := idxj r Y — D. The vector field D\ is said 
to be the total (or formal) derivative] the forms io 3 a are said to be contact (or 
Cartan) forms. Contact forms annihilate all vectors generated by D\. These 
are tangent to submanifolds of the form j r s(X) c J r Y. We will use iterated 
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total derivatives. Namely, if / : J r Y — ► M is a function, then we set At, a/ 
■=D x D„f. 

We have the following natural fibered splitting 

J r+1 Y T*J r Y = (j r+1 Y x T*X^j 8 C*[F], (1) 

where C*[Y"] :=imti;* is a subbundle of J r+ iY Xj r y T*J r Y, and is naturally 
isomorphic to J r +\Y x JrY V*J r Y (see pSlE^!)- 

The above splitting induces splittings in the spaces of forms [39] : here and in 
the sequel we implicitly use identifications between spaces of forms and spaces 
of bundle morphisms which are standard in the calculus of variations (see, e.g. 
[H EDI EH)- Namely, let A p be the sheaf of p-forms on J r Y. We introduce 
the sheaves of horizontal forms 'H p . +1 ri *- e - > °f fibered morphisms over 7r£ +1 
and 7r r of the type a: J r+1 Y -> A p T*J r Y and (3: J r Y -> A P T*X, respectively. 
Finally, for s < r we introduce the sheaves of contact forms C P s , i.e. , of fibered 
morphisms over 7rJ of the type a: J r Y — > A P C* [Y]). 

The splitting yields naturally the sheaf splitting 

p 

t=0 

Pull-back yields the inclusion A£ C W( r+lr) - The effect of © on is the 
following 

p 

A?C@C£-*AA*. (3) 
t=o 

Here, A' :=/i(A^) for < p < n and h is defined to be the restriction to A p 
of the projection of the above splitting onto the non-trivial summand with the 
highest value of t. Moreover, C P ~ l is the space of contact forms with values in 
C*[Y] (which is a bundle over J r+ {Y) and coefficients on J r Y. We define also 
the map v:=id— h. 

In other words, if a is a form on J r Y , then its pull-back (TT^ +1 )*a can be 
split into a part containing top degree horizontal forms and a part containing 
more contact factors (see, e.g. , 123125]). 

In coordinates this is achieved by means of the substitutions 

d^d\ dl^ul+yl +x d\ 

which allow to express a in the basis (d x ,uj l cr ) at the cost of raising the order of 
jet. 

The above splitting induces also a decomposition of the exterior differential 
on Y, (tt^ +1 )* o d = dtf + dy, where du and dy are called the horizontal and 
vertical differential, respectively The action of dn and dy on functions and 
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1-forms on J r Y uniquely characterizes dn and dy (see, e.g. , [33] for more 
details). In particular, we have the coordinate expressions 

d H f = D x fd x = (d x f + yl +x d?f)d\ d v f = dffui, 
d H d x =0, d H dl = -<4 +A A d x , d H uj l rT = -ujl +x Ad x , 
d v d x = 0, d v d l a = dl +x A d x , d V Lul = 0. 

Any fibered isomorphism F: Y — > Y over idx admits a prolongation to a 
fibered isomorphism J r F: J r Y — ► J r Y such that J r F o j r s = j r (F o s). 

A vector field ( on Y is said to be vertical if it has values into VY. A 
vertical vector field can be conveniently prolonged to a vertical vector field 
j r t;: J r Y — > VJ r Y. The vector field j r £ is characterized by the fact that its 
flow is the natural prolongation of the flow of £. In coordinates, if £ = we 
have j r £ — Da-^ 1 df , < |<x| < r. Again, if the order of prolongation needed in 
formulae for £ is clear from context, then we simply denote the prolongation of 

Let a £ C' ® A". Then we can interpret a as the differential operator 

where x$ is the space of vertical vector fields £ : Y — > VY and J denotes the inner 
product. This is an operator in total derivatives and its coordinate expression 
is V a (£ l <9;) = Drr^af vx, where vx '■= dx 1 A ... A dx n is the local volume form 
induced by a chart (x x ) on X. We can form the adjoint operator V* : (A™)* (g) 
A™ = T r — > ® A™ . It admits an intrinsic definition (see 0] for details) . Its 
coordinate expression is V* (/) = (— l)^D rT {a^ /) to 1 ®vx- This notion can be 
extended to a form a S Cl®Cl® A": in this case we have V Q , V* : x$ — » Cq® A", 
with coordinate expression 

V£(f'0i) = (-l) W iM«S£'V®i'x- (4) 
2.2 Variational sequence 

We recall now shortly the theory of variational sequences on finite order jet 
spaces, as it was developed by Krupka in [2"3| . 

Denote by (dker/i) s the sheaf generated by the presheaf dker h. Set 0* := 
kerft, + (dkerh) s . Then the restriction of exterior differential yields a sheaf 
sequence — > 6* which is an exact subsequence of the de Rham sequence of 
differential forms on J r Y. |23) . Such a subsequence is made by forms that do 
not produce contribution in action-like functionals 15, 221 1261 07) . 

Definition 1 The quotient sequence 

o Ry _ A ° h Aj/ej h. ... £^zi Af/ef ^ Af+ 1 ^±l...Ao 

of the de Rham sequence with respect to the contact sequence is called the r-th 
order variational sequence associated with the fibered manifold Y — > X. Here 
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the integer N depends on the dimension of the Cartan distribution on J r Y |23| . 

□ 

The variational sequence is locally exact and, due to the 'abstract' de Rham 
Theorem, it has the same cohomology as the de Rham sequence (see |23p. 

Pull-back of forms through 7r^ +1 yields a natural inclusion of the r-th order 
variational sequence into the (r + l)-th order variational sequence. This allows 
us to represent any equivalence class [a] G Ajf /Ojf with a single morphism of 
bundles over jet spaces, at the cost of raising the order r. More precisely, the 
quotient sheaves / <d k in the variational sequence are represented as sheaves of 
bundle morphisms (defined on jets of order s > r) via the following intrinsic 
isomorphisms [BlAI^O) 

h : K h ;iQ k ; -> V r fc : [a] ^ h(a), k < n, (5a) 
I n+ i : A» +1 /e? +1 - V? +1 : [/?] » E m , (5b) 

i n+2 : A?+ 2 /e; i+2 - v;+ 2 : [ 7 ] -> h Ki) . (5c) 

Let us describe the above morphisms h(a), E^p), -^(7) an d spaces V£ , 
/l < n + 2. 

1. /i(a) is just the horizontalization of a; in the case k = n the morphism h(a) 
can be interpreted as a Lagrangian density [Tl 1X771 ITTTl f^l f^Tl [TTl PTH] . 

2. Eh(p) is the Euler-Lagrange morphism associated to h(f3). More precisely, 
it can be proved |19U22ll3*^] that any form h(/3) can be uniquely split into 
the sum 

(4 r + + 1 1 )*h(P)=E m -d H (p m ) (6) 

where E h{p) G C^ 2r 0) A V£ r and p h{0) G C ( 1 2r _ l r _ 1) A V£rJ v Here E h{0) 
and diiiphip)) are uniquely defined, but Ph(/3), the momentum, is not; see 
[T§| for a discussion. 

3. is the Helmholtz morphism 1211 BU] associated to the form h(^f). 
More precisely, it can be proved that any form h(^) can be uniquely split 
into the sum 

(7r r 2 ;+ 1 )*/ l (7)-^( 7 )-rf ff (gh( 7) ), (7) 

where H h(l) G C^+i.r) A C (V+i,o) A V 2r and q h(p) G C| r _! A Vg^ (here 
uniqueness is intended in the same way as above), with the additional 
condition that Hh(-y) is skew-adjoint in the first contact factor. 

We will use later on the special case: 

[(^r+Y)H = [dE hW ] G £„ +1 (A^+ 1 1 /e^ + + 1 1 ). In this case, if we set 
7] := £ti(/3), then Hd v can be introduced as the skew-symmetrization of the 
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morphism H dr] e CL.+i^r+i) ® C (4r+i,o) ® W 4r+i> which is characterized 
by 

£(i jS dv)=jSjH dri (8) 

(see [!n] for details). We recall that 7? and H have the same kernel; 
more precisely, Hdrj = if and only if Hd n = 21 . 

Let us recall the coordinate expressions. 

1. h(ct) = Avx, where A € C°°(J r +iY) is a 'special' polynomial in the 
derivatives of order r + 1 (in the sense of see also EH1) • 

2. Being locally = Bf ui l a A vx, we have the standard expression of the 
Euler-Lagrange morphism (see, e.g. , |19l 1221 135] ~) 

E hm = {-lYD„Bf u* Av x . 

3. In the simpler case [(t^ 1 )*^ = \dE, m ] G £„ + i(A^+ 1 1 /e^+ 1 1 ), that we 
will use later on, we have locally dEum = dfejU)^ A w J A ux, where 
ej = (-l^D^Bf, so that 

H d E h(P) = H% <4 <8> ^ J <8> fx , HdE h 0) = ujI A w j A w x , (9) 

|p|=Q V IPI y 

where (cr, p) denotes the union of the multi-indices <x and p (see [211 1391 
HO) : a local version has been also derived in |2*3ll55| '). 

Next, we interpret the above spaces , with h < n + 2. 

1. V* := AjS, fc < n. We recall that A* = /i(A^). So, V™ can be interpreted as 
the space of Lagrangians of order r + 1 which are polynomials of 'special 
type' [2Hj in the higher order derivatives; 

2. V™ +1 is the space of Euler-Lagrange morphisms associated to forms h{f3); 
it is a subspace of Cq A W^r+i 

3. V™ +2 is the space of Hclmholtz (or Hclmholtz-Sonin) morphisms associ- 
ated to forms h{^). 

We can read Ek through the above isomorphisms Ik- We obtain the exact 
sheaf sequence 

► M — A° A vl A . . . f=±i v; i+2 — . . . (10) 

It turns out [231221 that: 
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1. if /i G V*, with k < n — 1, and = /i(a) = /fc(a), with a £ Ajf, then 

= /i(<£a). Notice that /i« +1 *cZa) = h((d v + d H )(h(a) + v(a))) = 
Dff(h(a)), hence £^ is equal to dn up to a pull-back; 

2. if A € V™ then £ n (X) £ V™ +1 coincides with the standard higher order 
Euler-Lagrange morphism associated with the Lagrangian A. We will 
simply write £ instead of £ n ; 

3. if 1] € V" +1 then £ n+ x(rj) € V™ +2 is the Helmholtz morphism correspond- 
ing to the Euler-Lagrange morphism r\. The exactness of (|1(J|) implies that 
£ n +i(?7) = if and only if there exists (locally) a Lagrangian A € V" such 
that £(X) — T), i.e. , rj is locally variational. 

Definition 2 Let k < n + 1. We say elements /i 6 V* to be variational forms. 

Remark 1 We observe that the spaces Aj? / 0^: with fc > n + 2 do not have in 
the literature (to our knowledge) any interpretation in terms of standard objects 
of the calculus of variations. In any case, there is a representation Ik also for 
such quotient spaces 00] ■ It comes from the analogue representation for the 
variational sequence on infinite order jets 0] p. 192]. Furthermore, variational 
forms of degree k > n + 1 will not play any role in the rest of the paper. LI 

Remark 2 Due to £k ° £k-i — 0, Lagrangians of the form A e £„-i(V™ _1 ) 
are variationally trivial (i.e. , they have identically vanishing Euler-Lagrange 
equations). Our aim in this paper is to obtain an intrinsic model for the second 
and higher order variations. This is achieved in the literature in several ways, 
each of which differs from the others by a total divergence. Our viewpoint is 
different: we want to provide a model which does not suffer the above arbitrari- 
ness. 

□ 

In view of the above Remark, we factorize a part of the variational sequence 
as follows: 

f 

V™ V" +1 




where V™ := V™/£„_i(V™ _1 ), 7r is the quotient map and £ is the factor map. 

It is important to compute infinitesimal symmetries of objects in the varia- 
tional sequence. To do this, it is natural to look for vector fields X such that the 
standard Lie derivative operator Lx passes to the quotient into the variational 
sequence. 



M. Francaviglia, M. Palese, R. Vitolo 



10 



The jet prolongation jE of vertical vector fields S : Y — ► preserves the 
contact structure on jets. Hence, it is easy to see that the standard Lie derivative 
operator with respect to jE preserves the contact sequence too. This yields the 
new operator Cjs. on the elements of the variational sequence 

C jS : V^V^.r^I k ([L jS a}), (12) 

where [a] = Jr (t). The operator Cjs is said to be the variational Lie deriva- 
tive |13| . Such an operator allows us to recover several well-known formulae 
from the calculus of variations (see, e.g. , in a unique picture. We have the 
following expressions [T3| : 

1. ifO<p<n-land/x€ V*, then 

£j£l* = jt\dvW (13a) 

2. if p = n and A S V™, then 

£ jC A = £J£(A) + <Mj'£J?Wa); (13b) 

3. if p = n + 1 and r\ E V™ +1 , then 

Cm = £(S\ti) + j£\H*,. (13c) 

We remark that the operator Cjs: V" — * V™ factorizes to V™, since the stan- 
dard Lie derivative operator with respect to j'S commutes with . This fact, 
combined with l|13bjl . produces the further new quotient operator Cjs'- 

Cjs ■ V? — > V™ : A i— > £ jH A :=S J £ (A). (14) 

It is clear that this operator can be interpreted as the infinitesimal variation 
operator of a Lagrangian up to total divergencies. In the next section we provide 
a geometric model for such an operator and its iterated applications. 

3 Variations of forms 

We shall here introduce the variation of a form as infinitesimal multiparameter 
deformation. This is realized by taking iterated Lie derivatives of the form with 
respect to vertical vector fields. In this paper we shall consider only vertical 
variations. In fact, as it is easy to realize (see appendix 1), variations with 
respect to projectable vector fields do not change the results, since they just 
add an horizontal differential (a total divergence) which is not relevant for our 
results in view of Remark 

Let s : X — * Y be a section and i > 0. Let Si, . . . , Si be vertical vector 
fields on Y. Denote by ip l tk , with 1 < I < i, the flows generated by S;. Then 
the map 



r(ti,...,ti) ^^o.-.o^os 



(15) 
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is said to be the i-th variation of s generated by (Si, . . . , Sj). 

Let a € Aj; and let T be an z-th variation of the section s. Then the map 

A ' I[r; S](a) := dt 1 d *.dt i \tu...M=o(jrT(ti, . . . , Ufa) (16) 

is said to be the i-th variation of the form a along the section s. 

The following Lemma states the relation between the i-th variation of a form 
and its iterated Lie derivative. 

k 

Lemma 1 Let a : J r Y — > /\T*J r Y. Let F be an i-th variation of the section s 
generated by variation fields Si, . . . , Sj. Then we have 

A l [T;s](a) = (j r s)*L jrSl . . . L jr3i a. (17) 
Proof. By the above definitions, we have 

A4[r; s]{a) = dh dz ..dtM -^=° [UrWu ° ■ ■ • ° ^ ° s)ra] 
= {jrS) * dh. . .dt) ^---t*=^ jr ^ r ° • • • 

= (Jrs)*Lj r3l . . .Lj rS .a, 

where 4>t k are the vertical flows generated by S^ and we used the definition of 
prolongation of a vertical vector field (see section [5J. \qed\ 

From the above considerations it follows that the definition of variation can 
be given in terms of Lie derivatives with respect to prolongations of vertical 
vector fields, without any reference to a given section. 

Definition 3 Let S lv . . ,Sj be vertical vector fields on Y. Then the variation 
of a form a € Aj? is defined to be the operator 

A 1 [Si, . . .,Ei](a) :=L JrSl (L JrS2 (. . . (L jrSi a) ...)). 

Remark 3 The variation of an (n + /i)-form along a section s is clearly zero. In 
fact, the above Lemma shows that the variation is obtained through a pull-back 
on X via s and any (n + /i)-form on X is zero. So, the definition of variation 
is trivial for (n + /i)-forms. 

Nonetheless, we can use the more general Definition [3] of variation in all 
situations. But we loose the classical interpretation of variation as 'derivative 
along a parametrized family of sections'. D 

4 Variations of Lagrangians 

In this section we restrict our attention to variations of elements of the vari- 
ational sequence, and, in particular, to Lagrangians. Our task is to compute 
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variations 'up to variationally trivial forms'. In other words, we want to com- 
pute the quotient variations of Lagrangians in the variational sequence. Indeed, 
this is straightforwardly permitted by our definitions: i-th variations are made 
by Lie derivatives with respect to prolonged vertical vector fields, and they pass 
to the quotient in the variational sequence |13| . 

Moreover, we will devote special attention to variations along sections which 
are critical with respect to a given Lagrangian A, i.e. , sections s such that 
£(A) o js = 0. 

4.1 Quotient variation in the variational sequence 

Let a € Aj5, with k < n + 1, and Hi,. . . ,Hj be vertical vector fields on Y. We 
have 

MtA^Hi, . . . , H,]a]) = 4([£ iHl (■ ■ • (L jSi a) ...)]) 

= £ iHl (...(W fe (M))...), 

where C stands for the variational Lie derivative (see l|13al) . I)13b[) . l|13c|l L 

Definition 4 The operator 

tf*[S x , . . . ,Ei]I k ([a]) := A- Hl (. • • (A-e 4 4(H)) . . . ) 

is said to be the quotient variation of the variational form /^([a]) € with 
respect to the vertical vector fields Hi,. . . ,3j on Y". 
If s : X — > Y" is a (local) section of tt, then the form 

(js)*(S i [Z 1 ,...,Z i ]I k ([a\)) 

is said the quotient variation of /&([«]) along s. 

Remark 4 The above definition of quotient variation can be applied to all 
quotient spaces in the variational sequence (see Remark^. But, if we want to 
calculate variations of forms along sections, only the variations of forms [a] € 
with k < n are non-trivial (see Remark EJ). In this work we will just devote 
ourselves to variations of Lagrangians; nonetheless it would be interesting to 
investigate variations of elements in with k < n. 

4.2 Second variation 

In this subsection we fix a Lagrangian A S V™ and compute its second quotient 
variation along a critical section s. 

Let H: Y — > VY be a vertical vector field. It is natural to introduce an 
improved quotient variation on the space V™ (see l|ll|l). In fact, the operator 
CjS (|14f> is equal to the operator Cjs (|13b(l up to 'total divergencies', i.e. , up 
to £„_i-exact (variationally trivial) Lagrangians. We recall (see Eq. JTDJ) that 
E n -\ is equal to dn up to a pull-back, or up to higher order variationally trivial 
terms. 
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Definition 5 The operator 

^[S 1 ,... ! S i ][A]:=£ iHl (...(4'E,A)...) 

= Si J £(3 2 J £ (. . . Si-ijf (Si J S{\)) ...)) 

is said to be the quotient variation of the Lagrangian A. 
If s : X — > Y" is a (local) section of tt, then the form 

C7*)*(5 i |H 1 ,...,S i ][A]) 

is called the quotient variation of [A] along s. D 

Of course, the first quotient variation of A is just <r-[S][A] = 3 J £(X). 
Let Si , 3 2 be two vertical vector fields and let us consider the second quotient 
variation of A 

5 2 [3 X ,3 2 ][A] =3iJ£(S 2 j£(A)). 

Definition 6 We define the fibered morphism Sj J £(S 2 J £(X)) to be the Hes- 
sian morphism associated with the Lagrangian A. 

We state our main result: 

Theorem 1 The second quotient variation of a Lagrangian A along a critical 
section s is equal to either one of the following intrinsic bundle morphisms, 
which are self-adjoint: 

1. the differential V£(X) of £(X) along the fibres of tt + , also known as 
vertical differential: 

V£(X) : J 2r +iY -» V*J 2r +iY ® V*Y <g> A n T*X; 

2. the adjoint V£{X)* of the vertical differential: 

V£(\y : J 2r +iY -► V*J ir+2 Y ® V*Y ® A"T*X 

(see QJ. 

Proof. In coordinates we have 3j = 3^<9i, 3 2 = SjSj and A = Luj. Then 

^[s^sjw = (-i)Msi Ar(dns|£(^)i))u 

= (-l)l CT l5{ £ CT (d, CT 5 2 £ (£w),) w + (-l)Hs| Dri-i, dJ£{Lu)i) u>. 



(18) 

If s is a critical section, then the first summand of 5 2 [5i,5 2 ][A] in the right- 
hand side of l|18(l vanishes identically. The second summand admits an intrinsic 
interpretation. In fact, we have £ 2 {X) = = H £ ( X ) due to the property £ n +\ ° 
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£ n = of the variational sequence and the fact that H and H have the same 
kernel Hence we have, from the expression JSJ 

0=jH 1 JiS 2 J# £(A ) (19) 
= Hi D^-i dj£{\) t to - (-1)1^ D^H* dJ£(Lu)i) u, (20) 

so that the second summand of 1)18(1 is equal to the first summand of i)20[l . 
which is the vertical differential T^£(A) of £ (A) (also known as linearization 0]) 
contracted with the prolonged fields Si, 5%, namely 

Si J j~2 J V£(X) = Si D ff S J 2 dj£{\)i to . (21) 

A comparison of the coordinate expressions shows that the second summand of 
1)20(1 is equal to the adjoint of V£(X) (see (0J). More precisely, 

SiJjS 2 jV£(A) =Si AjE 2 jV£{\)*. SMI 



Corollary 1 The morphism V£(X) is symmetric along any critical section s, 
i.e. 

J 2 ,.+is*(Si J jS 2 J V£(X)) = j2r+is*(E 2 J jEi J V£(X)). 

Definition 7 We define the fibered morphisms: V£(X)* = V£ (A) to be the 

Jacobi morphism associated with the Lagrangian A. 

It is not difficult to check that the above definition recovers the definitions 
given by several authors up to 'total divergencies' (see, e.g. , |14l 1261 
Our formulation however has the following advantages: 

1. it is manifestly intrinsic (or covariant); 

2. it holds for Lagrangians of arbitrary order, while in literature only the 
first order case is usually treated; 

3. it allows an easy generalization to iterated variations of any order; 

4. it allows an easy generalization to all spaces in the variational sequence. 

Remark 5 We can compare our approach with the one of Crampin ^U]. In 
that paper the Poincare-Cartan form is used to achieve an intrinsic formula 
for the second variation of a first-order Lagrangian (7' = 1) in the case of one 
independent variable (n = 1). However, as is well-known, Poincare-Cartan 
forms are no longer unique in the case of many independent variables (n > 1) 
and higher order Lagrangians (r > 1), so that our approach seems to be more 
suitable in the general case. 



M. Francaviglia, M. Palese, R. Vitolo 



15 



5 Examples 

Here we show by simple but relevant examples that our definition of Jacobi 
morphism coincides with the standard one and it can also be applied to higher 
order Lagrangians. 

Example 1 (Metric Lagrangian). Here we shall derive the classical Jacobi 
equation for geodesies within our framework. Let (Q,g) be a Riemannian 
manifold, with metric tensor g = g ab d a ®d h . The Lagrangian for geodesies 
is A — \g ab {q)q a q b dt and the Euler-Lagrange equation is given by 

£(A)„ = -[W +r abc q b q c ] = , 

where T abc are Christoffel symbols of the first kind. 

The Jacobi equation is then obtained by evaluating the local coordinate 
expression for the Jacobi morphism given by (|21|l for £(X) a . It is easy to see 
that the Jacobi morphism for geodesies is in fact given in local coordinates by 

V£(X)* = -[d a g bc q b + d a T bdc q b q d }^l + 
+ [d a9b A b - d c9ab q b ]~ a 2 ^~l + 
d d 2 

a b g a cq^2 dr i s„c- ad 2 t -i- 

Taking into account the Euler-Lagrange equation we get finally the Jacobi equa- 
tion 

d d 2 

[d a gbc^ b de ~ d c T dea ]q d q e E'i - 2T cba q b —E c 1 - g ac -^~\ = , 

which can be recasted in the standard form 

V? Si + iftem(Si,7, 7) = , 

where 7 is any geodesic curve, denotes the second order covariant derivative 
along the curve 7 and i?iem(Si, 7, 7) is the Riemannian curvature tensor. This 
agrees of course with |S]. 

Example 2 (Hilbcrt-Einstein Lagrangian). Let dimX = 4 and X be ori- 
entable. Let Lor(X) be the bundle of Lorentzian metrics on X (provided 
that it has global sections). Local fibered coordinates on J2(Lor(X)) are 

1 g^v> g^v.ai giiv,op)- 

The Hilbert-Einstein Lagrangian is the form Xhe £ defined by Xhe = 
Lheu, were L^e = r \fg~- Here r : J2(Lor(X)) — * M is the function such that, 
for any Lorentz metric g, we have r o j 2 g = R, being R the scalar curvature 
associated with g and g the determinant of g. 

A direct computation of the Euler-Lagrange morphism shows that Ed\ HE = 
G := Ric — I Rg 6 Ch Q ^ A 7i|, Ric being the Ricci tensor of the metric g. 
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The Jacobi equations for the Hilbert-Einstein Lagrangian can be then char- 
acterized as the kernel of the adjoint of the linearization morphism V£(Xhe)' 

+^r pA r^+(Gf+i s ^)r+v^v A r A +v ct ( 5 ^v A r A )- ( 5fv A ( 5 A ''v p rf;)] = o, 

which coincide with the classical variation of the Einstein tensor (see, e.g. , [H]). 

It is easy to realize that, along critical sections (i.e. solutions of the Einstein 
equations), the Jacobi morphism is in fact self-adjoint. This is also in accordance 
with [211. 

6 Conclusions 

We provided an intrinsic formalization of higher variations of a Lagrangian. In 
the case of the second variation, we gave a new interpretation of the Hessian 
and Jacobi morphism for Lagrangians of arbitrary order. 

Some problems remain open at this point and will be investigated in the 
future: 

• It would be interesting to derive a formula for higher order variations, as 
well as higher order analogues of Hessian and Jacobi morphisms. 

• It would also be worth to compute variations of all variational forms (not 
only Lagrangians). 

• There are branches of Quantum Field Theory in which higher order vari- 
ations play important roles, like e.g. the Batalin-Vilkoviski theory |17| . 
Such approaches still need a complete mathematical understanding. The 
above framework could be well-suited for that purpose: in |17| the Batalin- 
Vilkoviski theory is formalized through jet bundles. Moreover, higher or- 
der variations play in any case a role in the path integral approach to 
quantization. 

Acknowledgments 

Thanks are due to I. Kolaf and A. M. Verbovetsky for many useful discussions. 
The second author (M. P.) wishes to thank I. Kolaf for the invitation and the 
kind hospitality at the Department of Algebra and Geometry of the University 
of Brno in December 2000. 

Commutative diagrams have been drawn by Paul Taylor's diagrams macro 
package. 



M. Francaviglia, M. Palese, R. Vitolo 



17 



Appendix 1 

In this paper we considered only vertical variations. More general variations 
could be considered: in some problems of field theory (like the computation of 
conserved currents) it is interesting to consider Lie derivatives of Lagrangians 
with respect to protectable vector fields. A projectable vector field on Y is a 
pair (S, 5) such that H is a vector field on Y, S is a vector field on X and 5 is 
a bundle morphism over H. In coordinates, 5 = S a 9a + S l <9i, where S = S a 9a- 

However, such general variations do not modify the conclusions of our paper 
in a significant way: their contribution to the variation of a Lagrangian is in fact 
a total divergence, so that it has obviously to be neglected in our scheme. This 
fact is, of course, well-known (see, e.g. , ^H])? but we recall it here for the sake 
of completeness. The equation (|13b(l takes the more general form Cj r s,{\) = 
SyJ£(A) + d H (jr^v\Pd v \ +SJA) (see, e.g. , 13 ), where Zy :=w(3): J{Y -> 
VY is the vertical part of S. In coordinates, Sy = (S* — y\E x )di. 

It follows that our results of section 01 hold practically unchanged by just 
replacing vertical vector fields with vertical parts of projectable vector fields. 

Appendix 2 

Here we shortly recall the formulation of variational problems on jet spaces 
[II3imiI3E2ElE3EI3EE!to help the reader to connect the purely differ- 
ential setting of variational sequences with the classical integral presentation. 
Nonetheless, we stress the two approaches (differential and integral) are com- 
pletely independent, even if the latter provided the motivation to the former 
from an historical viewpoint. 

Suppose that an r-th Lagrangian A € is given. Then the action of A on 
a section s : U — > Y (U is an oriented open subset of X with compact closure 
and regular boundary) is defined to be the real number 

/ Urs)*X. 
Ju 

A vertical vector field £ : Y — > VY defined on ir^ 1 (U) and vanishing on 7r _1 (9E/)| 
is said to be a variation field. A section s : U — ► Y is said to be critical if, for 
each variation field with flow (fit, we have S JuiJr'fit ° jrs)*X — 0, where S is the 
derivative with respect to the parameter t and J r <fit : J r Y — > J r Y is the jet pro- 
longation of the flow (fit . It is easy to see that the previous integral expression is 
equal to J u (j r s)*L Ur X = for each variation field u, where u r : J r Y — > VJ r Y 
is the r-th jet prolongation of u (see the first section). Using equation i!3bl 
together with L Ur X = i Ur d\ and Stokes' Theorem, we find that the above equa- 
tion is equivalent to J u (j 2r s)* '{i u Ed\) = for each variation field u. Finally, by 
virtue of the fundamental Lemma of calculus of variations the above condition 
is equivalent to {j2rs)* Ed\ — 0, or, that is the same, Ed\ o j2 r s = 0. 

Now the reason of the choice of the sheaf 9^ (for < k < n) as the first 
non-trivial sheaf of the contact subsequence is clear: for k — n the sheaf 0™ is 
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made by forms which do not contribute to the action. As for the sheaf it 
is easily seen that this is precisely the sheaf of forms that give no contribution 
to the integral f u (j2 r s)*i u Ed\ when added to E d \. 
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